Properties of hadron screening masses at small baryonic density 
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The properties of hadron screening masses around the deconfinement phase transition at finite baryonic density 
can be studied by evaluating the Taylor coefficients with respect to the iso-scalar (ps = A*« = fid) and iso- vector 
(fiv ~ flu = — fid) chemical potentials, where fi u and fi^ are u and d quark chemical potentials. We simulate 
2-flavour lattice QCD with staggered fermions on a 12 2 x 24 x 6 lattice with ma = 0.05 and 0.10 and report 
investigations of nucleon, pseudo-scalar and vector mesons. We present new, strong evidence that in the confining 
phase, the screening masses at [i — have significant temperature dependence, but the effect of fis is very small. 
Above the critical temperature, the second derivative terms of mesons rapidly increase as contrasted to the case 
of baryon. We also study the responses of the screening masses to an iso-vector chemical potential and discuss 
some of the issues related to the properties of hadron masses at finite fx. 



1. Introduction 

There have been extensive studies of properties 
of hadron masses at finite temperature using phe- 
nomcnological models UI2| and lattice QCD sim- 
ulations |3I4I5| . while little is known about their 
behavior in the medium. Although many model 
studies and simulations of two-color QCD pro- 
duced interesting outcomes 6 , no lattice QCD 
study of hadron pole masses at finite baryon 
density has yet been performed. This is due 
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to the difficulty in performing numerical simula- 
tions with a complex fermion determinant. See 
Refs. [70. 

The chemical potentials of u and d flavors can 
be combined into iso-scalar and iso-vector (also 
called isospin) chemical potentials as fs = fu = 
fid and fly = fi u = —fid, respectively. In this 
letter we study the response of hadron screen- 
ing masses with respect to both fis and fly- I n 
addition to the pseudo-scalar (ir + (ud)) and vec- 
tor mesons, we also study the chemical potential 
response of the nucleon mass (n(udrf)). To our 
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knowledge, this is the first attempt to investigate 
baryon masses at finite baryonic density in SU(3) 
lattice QCD. 

Our analysis focuses on the behavior in the 
vicinity of zero chemical potential = 0. This is a 
region which currently running relativistic heavy- 
ion collision experiments at RHIC explore. In ad- 
dition to the usual, baryonic (iso-scalar) chemical 
potential [as, we also present data for the iso- 
vector chemical potential [iy case. Incident nu- 
clei in high energy heavy ion collisions are usually 
neutron rich, and therefore the [iy dependence 
reported here may be observed in these experi- 
ments. 

Using a Taylor expansion method at fixed tem- 
perature T, coupling f3 and bare quark masses, 
the screening mass is expanded at // = in the 
form, 
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This allows to perform numerical simulations 
with standard methods: we determine the Tay- 
lor coefficients at \i = 0, by measuring hadron 
propagators and their derivatives in [i = lattice 
QCD calculations |1U|. 

2. Lattice setup 

We study effects of the chemical potential on 
hadron screening masses with staggered fermions. 
The numerical simulations are carried out on a 
lattice of size 12 2 x 24 x 6, at eight values of (3 
from 5.30 to 5.65, with quark masses ma — 0.05 
and 0.10. Preliminary results were presented in 
Ref. jlOlllj . Earlier investigations were limited 
to smaller lattices 16 x 8 2 x 4 and small num- 
ber of j3 values. We report here new calculations 
around the phase transition, using lattices with 
a larger spatial volume and temporal extension 
Nt = 6, which provides us with data nearer to 
the continuum limit. 

The R-algorithm was used to generate 1000 
configurations each separated by 10 unit-length 
trajectories, with the molecular dynamics step 



size set to 0.02. Anti-periodic (periodic) bound- 
ary conditions were imposed in temporal (spatial) 
directions on the fermion fields. 

We deduced the scale from previous estimates 
of lattice spacings a for various values of (3 and 
quark masses determined by Tamhankar and Got- 
tlieb 12 . In the present analysis the lattice spac- 
ing a varies between 0.09 and 0.26 fm. 

The critical value of the coupling (/J c ) is esti- 
mated from the peak of the Polyakov line sus- 
ceptibility, giving (] c — 5.499 for ma — 0.05 and 
j3 c = 5.568 for ma = 0.10. The values of (3 studied 
here are listed in Table 1 together with a rough 
estimate of T/T c . 

To improve the ground state signal, we intro- 
duced a "corner" wall source. This source is not 
gauge invariant, so we fixed the link variables in 
the spatial slice z to the "Coulomb" gauge for 
each configuration. 

The Taylor coefficients in Eq. are obtained 
by taking derivatives of hadron correlators, for 
which we need the traces of various fermionic op- 
erators. To evaluate these traces, we used 200 
"noise" vectors of complex random numbers. 

3. Hadron screening masses 

Figure shows the temperature dependence of 
hadron screening masses obtained from fitting the 
correlators to a double-exponential form, 2 

C p (z) = C£ ( e - A i*+ e -*i (2) 

c N {z) = c'l (e~< i + {-iy e -™" 

+ C£ ((-1)* e-*i' a + e - *? t^"*)) , 

where N z is the length of the lattice in the z- 
direction. We write the correlators in terms of 
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dimensionless lattice variables, by scaling mi, m 2 
and z according to their canonical dimensions. 

We start our analysis with the hadron screen- 
ing mass which is obtained fitting parame- 
ter in each channel. Its value is the zeroth order 
term in the Taylor expansion of the hadron mass, 
i.e., M(fi = 0). For the pseudo-scalar meson, we 
do not consider a contribution with the alternat- 
ing sign. In this case, the two-particle fit includes 
an excited state with the same quantum numbers. 
For the vector meson and nucleon the two-particle 
fit includes one particle of each parity We 
have observed that a single-exponential fit is not 
good enough to extract reliable results in the de- 
confinement phase, where parity partners are ex- 
pected to be degenerate. 

The \x = screening masses divided by T are 
shown in Fig. ^ As expected, the pseudo-scalar 
and vector mesons become degenerate above T c 
because of chiral symmetry restoration. With 
increasing temperature, the screening masses of 
mesons (baryon) approach 2irT (3itT) which cor- 
responds to free quarks having the minimal mo- 
mentum 7rT 14 . The vector meson screening 
mass approaches this asymptotic value at a sur- 
prisingly low temperature. 
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Figure 1. The screening masses of pseudo-scalar 
(PS), vector (V) mesons and nucleon (N), divided 
by the temperature T. The dashed (dot-dashed) 
line represents the free-quark limit 2tt (Sir). Open 
symbols correspond to the quark mass ma — 0.05 
and filled ones to ma — 0.10. 



Table 1 

Values of [3 studied here and corresponding tem- 
peratures. 





ma = 0.05 


ma = 0.10 


a 


T/T c 


T/T c 


5.30 


0.58 


0.51 


5.35 


0.66 


0.56 


5.40 


0.76 


0.64 


5.45 


0.87 


0.74 


5.50 


1.01 


0.85 


5.55 


1.15 


0.95 


5.60 


1.31 


1.08 


5.65 


1.54 


1.24 



3.1. The pseudo-scalar meson 

The first derivatives of both the pseudo-scalar 
and the vector meson correlators with respect to 



the chemical potential vanish. In the case of the 
iso-scalar chemical potential, this can be easily 
seen from the symmetry fig — > — fig of the corre- 
lator. See Ref. |10lll| . The first derivatives with 
respect to the iso-vector chemical potential are 
numerically also consistent with zero. 

The second order responses of the pseudo- 
scalar meson correlator are shown in Fig. [21 In 
the confined phase, the response to fis is small: 
the pion is protected from acquiring a large mass 
because chiral symmetry is broken. In the de- 
confinement phase, the response increases with 
temperature and therefore the pseudo-scalar me- 
son becomes heavier when fig is finite. In the 
same figure we show data for two different quark 
masses ma = 0.10 and ma — 0.05. With de- 
creasing quark mass, the n$ response increases 
for T > T c . 

The response with respect to \xy shows a very 
different behavior. In the confinement phase it is 
negative, while above T c it approaches zero. If 
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we insert the values into Eq. QJ, neglecting the 
0[(/i/T) 3 ] terms, the mass becomes zero around 
fiy ~ M v , which is the expected critical value 
of the iso-vector chemical potential where a pion 
condensate forms .15 . This indicates that the 
pion mass dependence on fiy is rather well de- 
scribed by our truncated quadratic Taylor expan- 
sion in the whole confining phase. Equivalently, 
the second-order response can be compared with 
-2T/M PS (Fig.EJ), i.e. the value it would take if 
the quadratic ansatz were exact. 
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indicates that the vector meson screening mass 
above T c increases in the medium, assuming it 
still represents relevant degrees of freedom there. 

The second derivative with respect to fiy below 
T c shows similar behavior as the pseudo-scalar 
case. Above T c it approaches zero quickly and 
has little quark mass dependence. 
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Figure 3. The second order response of the vector 
meson. Open symbols correspond to the quark 
mass ma = 0.05 and filled ones to ma = 0.10. 



Figure 2. The second order response of pseudo- 
scalar meson. The open symbols correspond to 
the smaller quark mass ma = 0.05, the filled ones 
to ma = 0.10. The dotted (dashed) line stands 
for —2T/Mps, where Mps is the pseudo-scalar 
meson mass at ma = 0.05 (ma — 0.10) obtained 
from Fig.l. 



3.2. The vector meson 

As shown in Fig. |3 the situation with the 
vector meson is quite similar to the case of the 
pseudo-scalar meson described above. In the con- 
finement phase, the response to fj,s is almost zero, 
and increases in the deconfincmcnt region. This 



3.3. Baryon (nucleon udd) 

Unlike the case of mesons, the first order re- 
sponse of the nucleon does not vanish for either 
Us or \iy (see Fig. 0}. 

The derivative with respect to the iso-scalar 
chemical potential is positive and close to zero 
at T < T c . It increases quickly in the deconfinc- 
ment phase, and more so for smaller quark mass. 
The response to [iy is again negative. 

As for the second order response (Fig. |SJ, it is 
positive under the influence of fiy and negative 
in the case of /ig. The difference between the two 
quark mass cases is very small till 1.1T C . 

The second order term in the Taylor expansion 
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is of the same order as the first one, and gives a 
contribution of opposite sign. Consequently, the 
nucleon mass is less sensitive to chemical poten- 
tial than meson masses. 

4. Concluding Remarks 

We have analyzed the behavior of the first and 
second order derivatives of the meson and baryon 
screening masses with respect to the chemical 
potential for two quark masses ma = 0.05 and 
ma = 0.10. 
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Figure 4. First order response of the nucleon 
(neutron). Open (filled) symbols correspond to 
the quark mass ma = 0.05 (ma = 0.10). 



Our study shows that in the vicinity of zero iso- 
scalar chemical potential, the hadron screening 
mass increases as function of both temperature 
and density. At high temperature, the second or- 
der derivative with respect to fis increases and 
that w.r.t. [iy approaches zero for both pseudo- 
scalar and vector meson cases. At very high tem- 
perature the meson propagators should be super- 
position of free quarks. If we introduce an imag- 
inary chemical potential in momentum space 



T3 



3 
H 



1.5 
1 

0.5 


-0.5 



~~ 1 1 1 r 

H=Mv 



).4 0.6 0.8 1 1.2 1.4 1.6 

T/Tc 



Figure 5. Second order response of the nucleon 
(neutron). Open (filled) symbols correspond to 
the quark mass ma = 0.05 (ma = 0.10). 



they are written schematically as 

X] (7o(sin(p[ ) 1) + fj, Iu ) + fsinp^y 

p(l),p(2) 

r^osin^+^+fsinp^) V, (3) 

where po = rnrT with n = ±1,±2, The 
baryon propagator is constructed in the same 
manner. 

Such a "meson" propagator does not of course 
have a pole, but its behavior gives us a hint to 
understand what happens at finite chemical po- 
tential. In case of the iso-scalar chemical po- 
tential, the lowest contribution to the damping 
in the z-direction comes from | — ttT + m\ + 
| — (— 7rT + ~ 2(t:T — hi) if /!/ is posi- 
tive and smaller than 2-kT. For the iso-vector 
case, we have \irT — + \ — ttT — ~ 2nT. 
Therefore when T increases above T c , the meson 
mass reaches 2ttT showing little dependence on 
the imaginary iso-vector chemical potential, while 
with the iso-scalar chemical potential, the mass 
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decreases and may behave as M(0) — Cpj with 
C positive. If we now go back to a real chemical 
potential, M(0) + Cp 2 , i.e., the second derivative 
with respect to ps is positive, while that w.r.t. 
pv approaches zero as seen in Figs |21 and |3| 

The second response of the vector meson with 
respect to ps above the phase transition tempera- 
ture is similar to that of the pseudo-scalar meson. 
This indicates that when the chiral symmetry is 
restored, not only the masses of n and p, but also 
their response to the chemical potential is similar. 

Although the quark mass dependence below T c 
is different, the gross features of the response to 
pv of the p meson are similar to that of the tt. 
This may indicate an interesting phenomena, i.e., 
in an experimental environment where the num- 
ber of u quarks is not equal to that of d quarks, 
the vector meson mass may shift. 

There have been many analyses of the fi- 
nite iso-vector chemical potential system by NJL 
model |18I19| , chiral lagrangians .15 , random ma- 
trix model [20| ) the strong coupling expansion [2~T] 
and lattice simulations |17I22I23I24I25| . Son and 
Stephanov first pointed out that at py — m, 
and T = 0, there is a phase transition into a 
pion condensation phase where 7r becomes mass- 
less According to these analyses, the lead- 
ing iso-vector chemical potential effect of the pion 
pole mass is a linear in py. Our simulations give 
an insight into the region of finite temperatures 
where we observe that the first response of the 
screening mass w.r.t. py vanishes. 

In Ref. ^Sj) QCD is studied at finite py in the 
low and high py limits. The calculation of the 
neutron mass is performed using baryon chiral 
perturbation theory. 

M = M(p = 0) - ^ cos(a), (4) 

where cos a = 1 for pv < m^ and cos a = m„/fj, v 
for pv > ra v . This may be consistent with the 
negative first derivative obtained here. 

The value of the first response of the nucleon 
with respect to the iso-scalar chemical potential 
(Fig. is very small but positive. If this is defi- 
nitely non-zero, it means that at low but positive 
baryon density environment, the nucleon screen- 
ing mass is heavier, i.e., the baryon exchange 



mechanism may be suppressed. 

The present simulation was performed with two 
fixed values of ma =0.10 and 0.05. In order to 
compare with experimental data, we should ex- 
trapolate them to the chiral limit ma — > 0. We 
plan to study the screening masses at finite T and 
p along a line of constant physics corresponding 
to the real (two flavor) world. 
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